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Spontaneous reconnection at a separator current layer. 

II. Nature of the waves and flows 

J. E. H. Stevenson,^ and C. E. Parnell^ 

Abstract. Sudden destabilisations of the magnetic field, such as those caused by spon¬ 
taneous reconnection, will produce waves and/or flows. Here, we investigate the nature 
of the plasma motions resulting from spontaneous reconnection at a 3D separator. In 
order to clearly see these perturbations, we start from a magnetohydrostatic equilibrium 
containing two oppositely-signed null points joined by a generic separator along which 
lies a twisted current layer. The nature of the magnetic reconnection initiated in this equi¬ 
librium as a result of an anomalous diffusivity is discussed in detail in Stevenson and 
Parnell [2015]. The resulting sudden loss of force balance inevitably generates waves that 
propagate away from the diffusion region carrying the dissipated current. In their wake 
a twisting stagnation-flow, in planes perpendicular to the separator, feeds flux back into 
the original diffusion site (the separator) in order to try to regain equilibrium. This flow 
drives a phase of slow weak impulsive-bursty reconnection that follows on after the ini¬ 
tial fast-reconnection phase. 


1. Introduction 

It has long been recognised that magnetic reconnection 
generates waves and flows since the magnetic energy released 
by reconnection not only leads to direct heating, but also ac¬ 
celerates populations of particles and the bulk plasma. In¬ 
deed, the simple order-of-magnitude estimates of the plasma 
behaviour in a steady two-dimensional (2D) magnetohy¬ 
drodynamic (MHD) reconnection scenario [Parker, 1957] 
revealed that reconnection outflows can be at the Alfven 
speed. Petschek [1964] recognised that such fast outflows 
could lead to shocks being created in the outflow regions, 
and developed a steady 2D MHD reconnection model incor¬ 
porating shocks producing both fast reconnection, as well 
as additional heating on top of that due simply to Ohmic 
dissipation alone. 

Many modifications have been made to these models with 
numerous more complex 2D reconnection scenarios proposed 
[see, for example, the reviews by Priest and Forbes, 2000; 
Biskamp, 2000]. In particular, with the ability to perform 
large-sale numerical experiments, 2D reconnection is now 
modelled using a wide range of approaches including MHD, 
Hall-MHD, multi-fluid, hybrid and particle-in-cell [e.g. Birn 
et al, 2001]. The addition of extra physics beyond MHD 
means that instead of just MHD waves (Alfven, fast mag- 
netoacoustic and slow magnetoacoustic) being found, there 
are also higher frequency waves such as Whistler waves 
[e.g., Drake et al, 1997; Fujimoto and Sydora, 2008] and 
ion/electron cyclotron waves [e.g., Hoshino et al, 1998; 
Arzner and Seholer, 2001], which can have a significant ef¬ 
fect on characteristics such as the onset time and the rate 
of reconnection. 

Reconnection occurs in many geophysical situations, e.g., 
solar flares, coronal mass ejections, substorms, interactions 
between planetary magnetospheres and interplanetary mag¬ 
netic fields (IMF). Observational evidence of fast flows from 
reconnection has been detailed for many years [Pasehmann 
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et al, 1979; Sonnerup et al, 1981; Gosling et al, 1986; 
Innes et al, 1997; Phan et al, 2000; 0ieroset et al, 2000; 
Yokoyama et al, 2001; Ko et al, 2003; Lin et al, 2005; Wang 
et al, 2007; Nishizuka et al, 2010; Milligan et al, 2010; Liu 
et al, 2010; Hara et al, 2011; Takasao et al, 2012; Savage 
et al, 2012; Watanahe et al, 2012; Cao et al, 2013]. Here, 
we are specifically interested in the nature of the waves and 
flows generated as a result of three-dimensional (3D) recon¬ 
nection: a topic which has not been widely studied. This 
may be due to several factors: (i) 3D reconnection has many 
differences to 2D reconnection and identifying exactly where 
the reconnection occurs and its nature are much harder to 
do in 3D than in 2D, also (ii) most models of reconnection 
are driven and so it is difficult to disentangle the waves gen¬ 
erated as a result of the reconnection from the flows driven 
by the boundary conditions. Due to these difficulties we 
constrain ourselves here simply to studying the waves and 
flows within a 3D MHD model generated as a result of recon¬ 
nection that occurs spontaneously in a magnetohydrostatic 
(MHS) equilibrium. 

In 2D there have only been a few models specifically de¬ 
signed to investigate the MHD waves generated by reconnec¬ 
tion [Longeope and Priest, 2007; Puentes-Fernandez et al, 
2012a, b; Longeope and Tarr, 2012]. Here, we briefly discuss 
the results of a 2D MHD model involving undriven (sponta¬ 
neous) reconnection occurring in a high-beta plasma, whose 
approach we follow to investigate the MHD waves generated 
due to 3D reconnection in a high-beta, MHD scenario. 

In order to study the nature of the MHD waves generated 
from 2D X-point reconnection, Puentes-Fernandez et al 
[2012a] used the approach of first forming a MHS equilib¬ 
rium with a current layer about a 2D X-point, before study¬ 
ing the reconnection and associated waves at the null em¬ 
bedded in a high-beta plasma. In Puentes-Fernandez et al 
[2012a], to trigger reconnection in the current layer, which 
may arise, for instance, as a result of micro-instabilities, an 
anomalous diffusivity was introduced. The addition of an 
anomalous diffusivity term, which acts only where the cur¬ 
rent is greater than a set value, leads to the current layer 
(and not the enhanced current along the separatrices) being 
diffused rapidly. 

Waves, launched from the diffusion site at the fast and 
slow magnetoacoustic speeds, travel outwards leaving a stag¬ 
nation flow pattern behind in their wake [Puentes-Fernandez 
et al, 2012a]. This flow is created because the system tries to 
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restore the equilibrium that has been lost, as a result of the 
reconnection, by rebuilding the current layer, but this sim¬ 
ply drives further reconnection. The magnetoacoustic waves 
carry current away from the current layer and propagate en¬ 
hancements/deficits of plasma pressure in the outflow/inflow 
regions. Since the fast and slow magnetoacoustic speeds are 
very similar in a high-beta plasma the outward propagating 
waves maintain an elliptical shape, although the major axis 
of the ellipse switches over time as the speed of the waves 
is quite different in the inflow and outflow regions. R was 
found that most of the reconnection in this high-beta case 
occurred during an initial rapid diffusion phase (which was 
followed by a second slow reconnection phase driven by the 
flows left in the wake of the waves). 

An identical experiment was run, but with a surround¬ 
ing low-beta plasma [Fuentes-Fernandez et al.^ 2012b], al¬ 
though, at the null itself and in its immediate vicinity the 
plasma is high-beta, as it must be by the definition of a 
null. In contrast to the high-beta case, there are distinct 
differences in the propagation of the magnetoacoustic wave 
pulses because the fast and slow speeds are distinct in a low- 
beta plasma. Additionally, in the low-beta case [Fuentes- 
Ferndndez et a/., 2012b] most of the reconnection occurred 
in the second phase, as opposed to the first, through an 
impulsive-bursty reconnection regime. Impulsive-bursty re¬ 
connection is not achievable in the high-beta experiment 
due to the low magnitudes of the forces left in the wake of 
the propagating waves. These differences in flows are high¬ 
lighted by the fact that the amplitude of the propagating 
waves of the high-beta case are 10^ times smaller than those 
of the low-beta case. 

As already briefly mentioned, reconnection in 3D is fun¬ 
damentally different to that in 2D and it can occur at topo¬ 
logical or geometrical features of a magnetic field [Hesse and 
Sehindler^ 1988; Sehindler et a/., 1988]. In this paper, we 
focus on the reconnection which occurs at a topological fea¬ 
ture called a separator, since such features have been shown 
to be prime locations for reconnection [e.g., Sonnerup, 1979; 
Lau and Finn, 1990; Longeope and Cowley, 1996; Galsgaard 
and Nordlund, 1997; Galsgaard et al., 2000; Longeope, 2001; 
Pontin and Craig, 2006; Priest et al., 2005; Haynes et al, 
2007; Parnell et al., 2008; Dorelli and Bhattaeharjee, 2008; 
Parnell et al., 2010; Komar et al., 2013; Stevenson et al., 
2015]. 

Only generic separators exist for more than an instant 
in dynamic 3D magnetic fields (separators formed by the 
intersection of the spines, or one spine and a separatrix sur¬ 
face, from two distinct 3D nulls are non-generic as any small 
perturbation in the field will destroy the intersection) and 
are formed by the intersection of the separatrix surfaces of a 
pair of 3D null points and so are special field lines that con¬ 
nect two null points [see Lau and Finn, 1990; Parnell et al, 
2010; Stevenson and Parnell, 2015, for a basic discussion on 
3D null points, separators and separator reconnection]. 

Although, it has been established that separators are a 
common topological feature found throughout the solar at¬ 
mosphere [e.g.. Close et al., 2005; Flatten et al., 2014; Par¬ 
nell et al., 2015], here we focus on those which may exist in 
the chromosphere, the Earth’s magnetosphere or other plan¬ 
etary atmospheres. Notably, in the dayside magnetopause, 
which is a high-beta plasma region [Trenehi et al., 2008], it is 
well known that reconnection will most likely occur along the 
line which separates the four topological flux domains con¬ 
taining the interplanetary magnetic field (IME), the closed 
terrestrial magnetospheric field lines and the open field lines 
that extend from the Earth out into the IME or extend down 
from the IME to the Earth. Erom a 2D perspective, these 
four flux domains only come together at a single point which 
must be a null point. In 3D, however, these four domains 
come together all the way along a line, the field line know 
as the separator, which crucially does not have zero mag¬ 
netic field all the way along it. The local magnetic field in 


planes perpendicular to a separator may be either X-type 
or 0-type in nature, as demonstrated both analytically and 
numerically by Parnell et al. [2010] and also found in Steven¬ 
son and Parnell [2015]. Thus, the name “X-line”, which has 
in the past been used to refer to such a line, is inappropriate 
and should simply be reserved for scenarios in 2.5D. 

On the dayside magnetopause many models have been 
formulated to predict where the location of this reconnec¬ 
tion occurs [e.g., Sonnerup, 1974; Gonzalez and Mozer, 1974; 
Alexeev et al, 1998; Moore et al, 2002; Trattner et al, 2007; 
Swisdak and Drake, 2007; Borovsky, 2008; Dorelli and Bhat¬ 
taeharjee, 2009; Borovsky, 2013; Hesse et al., 2013]. Komar 
et al. [2013] have mapped the dayside magnetopause separa¬ 
tors in global magnetospheric simulations for arbitrary clock 
angle of the IME and found that in all cases separators exist 
and are the locations at which the reconnection occurs. 

This paper is the second of a series. In the first paper 
[Stevenson and Parnell, 2015] the nature of the magnetic 
reconnection which occurred at a single-separator MHS equi¬ 
librium current layer, embedded in a high-beta plasma, was 
studied. In this paper, we focus on the properties of the 
MHD waves and flows generated as a consequence of sep¬ 
arator reconnection. In order to achieve this we study the 
local region about an isolated straight separator. Obviously, 
such an idealised scenario is unlikely to be realised in the 
solar chromosphere, or any planetary magnetosphere. How¬ 
ever, our qualitative results should be applicable in any of 
these scenarios (under the constraints of MHD) and our di¬ 
mensionless results may be scaled using dimensional factors 
to produce values that can be compared with results from 
larger-scale numerical models or observed quantities. 

We begin in Sect. 2 by briefly summarising the properties 
of the reconnection which are discussed in full in Stevenson 
and Parnell [2015] and then recap the details of the initial 
setup and the MHD code used to carry out the reconnec¬ 
tion experiment (Sect. 3). In Sect. 4 we analyse the waves 
launched, due to the reconnection, and then look at the 
transport of energy in the system (Sect. 5). Einally, we 
summarise our findings in Sect. 6. 

2. Nature of the Reconnection 

Stevenson and Parnell [2015] studied the properties of the 
reconnection which occurred at a separator current layer em¬ 
bedded in a high-beta plasma. They found that the recon¬ 
nection occurred in two distinct phase; a fast-reconnection 
phase (Ot/ < t < 0.09t/) in which 75% of the magnetic en¬ 
ergy was converted into internal and kinetic energy and a 
slow, impulsive-reconnection phase (0.09t/ < t < 0.76t/) in 
which only short-lived sporadic reconnection events occur. 
All times in the experiment discussed in Stevenson and Par¬ 
nell [2015], and discussed here, are normalised to the time 
it would take a fast-magnetoacoustic wave to travel along 
the MHS equilibrium separator from one null to the other 
{tf = 0.88). The experiment was stopped at t = 0.76t/ 
since the waves, launched from the diffusion site at the start 
of the reconnection experiment, neared the boundaries at 
this time. The sporadic reconnection events which occur 
in phase H were numerous enough such that the total flux 
reconnected continued to increase during this phase. The re¬ 
connection was observed to occur asymmetrically along the 
entire length of the separator and had a counter-rotating 
flow associated with it. 

Throughout both phases of the reconnection experiment. 
Ohmic heating dominated over viscous heating (and adia¬ 
batic cooling) due to the model being embedded in a high- 
beta plasma. The value of the plasma beta, which is high 
due to the presence of two null points in the system, in¬ 
hibits the waves from becoming large throughout the exper¬ 
iment. The strongest regions of reconnection occurred along 
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the separator away from the nulls and are associated with 
elliptical magnetic held lines in planes perpendicular to the 
separator. The length of the hrst (main reconnection) phase 
and the amount of magnetic energy which is converted was 
shown to be dependent on not only the size of the diffusivity, 
Tjd^ but also the size of the diffusion region, jcrit- 

The sudden loss of force balance in the MHS equilibrium 
caused by the anomalous diffusivity, leads to waves being 
launched from the diffusion site which travel outwards and 
cause changes to both the plasma and the magnetic field 
in their wake. In this paper, we detail the nature of these 
waves and fiows set up in the system by the reconnection, 
and analyse how they affect different plasma parameters. In 
the next section, we briefiy outline the properties of the MHS 
equilibrium current layer which we use as our initial state in 
our separator reconnection experiment and summarise the 
numerical model used for completeness. 

3. Initial MHS Equilibrium Current Layer 
and Numerical Model 

The initial state of our reconnection experiment is a MHS 
equilibrium which contains a twisted 3D current layer lying 
along the separator. Eig. la shows the MHS equilibrium 
skeleton along with the current layer (represented by an iso¬ 
surface of current drawn at jcrit = 10: the current above 
which the diffusivity is non-zero). This MHS equilibrium 
current layer was formed through the non-resistive MHD re¬ 
laxation of a non-force free magnetic field which contained 
two null points of opposite signs, whose separatrix surfaces 
intersected along the z-axis to form a generic separator. The 
full details of how a similar MHS equilibrium was formed 


and the properties of the separator current layer are given 
in Stevenson et al. [2015]. 



(c) -0.+ -0.2 0.0 0.2 0.4 


X 

Figure 1. Skeleton of the MHS equilibrium mag¬ 
netic field with (a) purple isosurface of j\\ = 10.0 and 
(b) yellow/blue isosurfaces of the pressure difference 
(p — Po, where po = 1.5) drawn at 70% of the maxi¬ 
mum positive/negative values. Also shown are the pos¬ 
itive/negative nulls (blue/red spheres) with associated 
spines (blue/red lines) and separatrix-surface field lines 
(pale-blue/pink lines) and the separator (green line, hid¬ 
den by the current layer in (a)). The solid pale-blue/pink 
lines indicate where the separatrix surfaces intersect the 
domain boundaries, (c) Perpendicular cut across the 
MHS equilibrium separator at z = 0.4 showing contours 
of the pressure difference with black and white lines show¬ 
ing contours of the magnitude of the current (|j|). Eight 
cyan asterisks are drawn in four positions on the edge of 
the white contour {jcrit = 10) which represents the edge 
of the diffusion region in this plane. The insert highlights 
the depth (d) and width (le) of the diffusion region in this 
plane. 

The plasma pressure in this equilibrium is such that pres¬ 
sure enhancements lie in cusp regions about the separator, 
and the pressure falls off away from here. Eig. lb shows the 
skeleton of the equilibrium field with yellow/blue isosurfaces 
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of the pressure difference (the pressure in the equilibrium 
state, p, minus the uniform pressure in the domain before 
the non-resistive relaxation takes place, po = 1.5). Con¬ 
tours of the equilibrium pressure difference are shown in a 
plane perpendicular to the separator at z = 0.4 in Eig. Ic. 
Over plotted here are black and white contours of the mag¬ 
nitude of the current |j| in this plane. The white contour is 
drawn at jcrit = 10 which is the value that represents the 
separator current layer in this plane. Eight cyan asterisks 
(shown in four positions) are plotted in Eig. lb which lie 
on the “edges” of the MHS equilibrium current layer on the 
current contour equal to jcrit = 10. The positions of these 
asterisks will be used (Sect. 4) to highlight the speed at 
which waves, launched by the onset of reconnection, move 
out from the separator current layer. 

This MHS equilibrium is used in Stevenson and Parnell 
[2015] as the initial state in a resistive MHD experiment us¬ 
ing the LareSd code [Arber et a/., 2001]. Line-tied boundary 
conditions are used (the normal components of the magnetic 
field, density and internal energy per unit mass attain min¬ 
ima or maxima on the boundaries) and the velocity is set 
to zero on the boundaries. Reconnection is triggered at the 
separator current layer, which existed in the MHS equilib¬ 
rium, through the use of an anomalous diffusivity which is 
zero unless the current is greater than the value, jcrit, where 
it takes the value rjd- As in Stevenson and Parnell [2015], 
we use jcrit = 10.0, such as to include the strong current in 
the separator current layer in the reconnection, but not the 
enhanced current on the separatrix surfaces, rfd — 0.001 and 
a constant background viscosity of = 0.01. 

Below, we detail the nature of the waves and flows which 
are created in the system as a result of the spontaneous 
reconnection discussed in Stevenson and Parnell [2015]. 

4. Propagation of MHD Waves 

Initially the plasma is in a MHS equilibrium, but as soon 
as the current in the current layer starts to dissipate, due to 
the onset of localised reconnection, waves are launched from 
the edges of the diffusion region (main current layer Eig. 1). 
These waves travel throughout the system communicating 
the collapse of the current layer and the resulting loss of 
force balance. In their wake, the magnetic field and plasma 
respond to these changes. In this section, we describe the 
nature of the waves launched and the resulting response of 
the plasma after they have passed. 

In order to investigate these waves we consider the per¬ 
turbed current, which we define as |j|-|jMifs|, rather than 
Cl'S |j —iMHs\ such that we can see both enhancements and 
deficits in the magnitude of the current (top row of Eigs. 2, 
3 and 4). We also examine the perturbed pressure, p-pMHS 
(middle row of Eigs. 2, 3 and 4) and the vorticity (bottom 
row of Eigs. 2, 3 and 4) with snapshots shown at three dif¬ 
ferent times to illustrate their behaviour. The three times 
which we show represent the experiment near the start of 
phase I (t = 0.019t/), at the end of phase I (t = 0.09t/) 
and about 75% of the way through phase H (t = 0.60t/). 
The waves launched are very small, with amplitudes of or¬ 
der 10~^ in current and of order 10~^ in pressure. The 
amplitudes of the waves are small in comparison to the 
MHS equilibrium values and reflect the size of the distur¬ 
bance that caused them, namely the dissipation of the cur¬ 
rent layer. As discussed in Stevenson and Parnell [2015], a 
high-beta plasma contains current layers that are fatter and 
have smaller maximum current than the identical low-beta 
plasma scenario, due to the stabilizing effect of the pressure 
gradient force in high-beta plasmas. See movies online for 
the full evolution of the current, pressure and vorticity in 
the planes shown in Eigs. 2, 3 and 4. 


Cuts in the ^ = 0.4 plane are plotted in Eig. 2 showing 
the behaviour of the perturbed current, perturbed plasma 
pressure and vorticity at three different times. In Eigs. 3 
and 4, the cuts are taken at the same three times, but are 
through the depth and across the width of the current layer, 
respectively. So, the vertical axis for these graphs runs along 
the z-axis from the lower to the upper null. However, the 
current layer is twisted and so the horizontal axis changes 
with distance along the separator, such that it always goes 
through the current layer depth (Eig. 3) or width (Eig. 4), 
as required. On all these graphs, we have plotted aster¬ 
isks, which at the start of the experiment are located in the 
z — A plane on the current contour |j| = jcrit, i-e., on the 
edge of the diffusion region (shown in the equilibrium field 
in Eig. Ic). We advect these asterisks through the depth 
and across the width of current layer in the z = 0.4 plane 
at the fast-magnetoacoustic speed, c/, both outwards away 
from, and inwards across, the current layer, in order to esti¬ 
mate the location, at a given time, of any fast-mode waves 
launched at the onset of the reconnection. Note, though, 
since the domain is three dimensional, the waves do not have 
to move in a planar manner, which may explain the small 
discrepancies between the wave fronts and the asterisks. 

4.1. Current and Pressure Perturbations 

The onset of reconnection causes a sudden deficit in cur¬ 
rent at the current layer, but the current does not disappear, 
instead it is carried away from the reconnection site by wave 
pulses that travel at the fast and slow magnetoacoustic wave 
speeds. Since the plasma is high beta the fast and slow mag¬ 
netoacoustic wave speeds are very similar in our experiment. 
Shortly after these waves have been launched (left-hand col¬ 
umn of Eigs. 2, 3 and 4 at t = 0.019t/), a distinctive wave 
pattern is setup in perturbed current and pressure. In the 
plane perpendicular to the separator (Eig. 2) this pattern 
is essentially the same as that found following the onset 
of reconnection at a 2D current layer [Puentes-Fernandez 
et a/., 2011]. The leading peaks of the current perturba¬ 
tions show enhancements in current followed by deficits, 
whichever direction they are travelling away from the re¬ 
connection site. However, the leading peaks of the pressure 
perturbations show an enhancement when travelling across 
the width (Eig. 4), but a deficit when travelling through the 
depth (Eig. 3) of the current layer. The former are launched 
from the narrow edges of the current layer, as if from a 
point, and propagate outwards, in a spherical-like manner, 
into the cusp regions, which in the MHS equilibrium, have 
larger pressure than the regions outside the cusps. 

Outside the cusps, the latter type of perturbations are 
launched (i.e., carrying a deficit in pressure) from the com¬ 
paratively wide edge of the current layer. These perturba¬ 
tions are more planar in nature, thus they move linearly 
outwards in a block from either side of the current layer. 
The left-hand columns of Eigs. 3 and 4 show that the per¬ 
turbations in current and pressure show a consistent be¬ 
haviour down the entire length of the separator current layer 
both through its depth and across its width. Remember 
though that the current layer is twisted and so these per¬ 
turbations emanate out along the separator forming an ex¬ 
panded helical pattern. This is shown in a movie, available 
online, where isosurfaces of the pressure difference (drawn 
at p — Po = 1 X 10~^ in yellow and p — po = —1 x 10~^ in 
blue) evolve through the reconnection experiment. In this 
movie, we are looking down on the 3D model and so only 
the negative null (red sphere) is seen but the spines of both 
nulls are visible. 



STEVENSON, J. E. H. AND PARNELL, C. E.: IL NATURE OE THE WAVES AND ELOWS 


X- 5 


0.0010 

0,0005 


0.0000 


-0.0005 


- 0.0010 


1x10"* 


5x10-® 


0 


-5x10"® 


-1x10“^ 


0.030 


0.015 


0.000 


- 0.015 


- 0,030 





Figure 2. Contours of |j| — HmhsI (first row), p — pMHS (second row) and Uz (third row) in the plane 
= 0.4 across the separator, at t = O.OlOt/ (first column), t = 0.09t/ (middle column) and t = O.GOt/ 
(last column). Asterisks, which initially lie on the edge of the diffusion region, as shown in Eig. Ic, 
move at the fast-magnetoacoustic speed, c/(x, y, z, t). Over plotted on the bottom row of graphs are 
arrows (normalised to the maximum value of the magnitude of the velocity in the domain at t = O.GOt/, 
|v| = G X 10“^) that display the direction of Vx and Vy in = 0.4 plane. As time increases, so do the 
dimensions of the planes. 
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Figure 3. As for Eig. 2, but instead showing the perturbations in a vertical surface that crosses the 
depth of the current layer at right angles to its width. Here, the arrows (normalised to the maximum 
value of the magnitude of the velocity in the domain at t = 0.60t/, |v| = 6 x 10~^) display the direction 
of Vr and Vz in this plane in the last row. The arrows are coloured black where d<0 and are coloured 
green where d>0. 
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Figure 4. As for Eig. 3, but instead showing the perturbations in a vertical surface that crosses the 
width of the current layer at right angles to its depth. Here, the arrows are coloured black where w<0 
and are coloured green where ie>0 and are normalised as in Eig. 3. As in the previous hgures, as time 
increases, so do the dimensions of the planes. 


As these perturbed pulses travel further from the re¬ 
connection site they are followed by a second set of 


pulses, which show the same basic behaviour. These 
pulses are the ones that were launched at the same 
time as the lead pulses, but travelled inwards across 
the current layer, rather than outwards from it. Nat- 
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urally, therefore, for the wave pulses travelling outwith 
the cusps, the following planar pulses are very close be¬ 
hind the lead planar ones since they simply cross the 
(thin) depth of the current layer (right-hand column of 
Figs. 2 and 3). In the cusps themselves, the following 
pulses that leave the narrow edges of the current layer 
have to cross the entire width of the current layer and, 
therefore, these are much further behind the leading 
point-like pulses (right-hand column of Figs. 2 and 4). 

In Fig. 5, time slices of the perturbed current and 
pressure through the depth and across the width of the 
current layer in the 2 ; = 0.4 plane are plotted. These 
show the wave pulses that travel from the edges of the 
diffusion region and match well with the over plotted 
lines that indicate the speed that fast-magnetoacoustic 
waves, launched both inwards and outwards from the 
edges of the diffusion region, would travel at. 




Figure 5. Time slices of the perturbed current (up¬ 
per row) and the perturbed pressure (lower row) plotted 
through the depth (first column) and across the width 
(second column) of the current layer, in the plane at 
2 : = 0.4. The black dashed lines highlight where the 
first phase ends and phase II begins. The green/white 
lines start on the edge of the current layer and represent 
a wave travelling at the fast-magnetoacoustic speed. 


4.2. Steady Flows 

In the wake of these perturbations, deficits of cur¬ 
rent exist both outside and within the cusps. Deficits 
of pressure exist only outwith the cusps, while enhance¬ 
ments of pressure inside the cusps form locally around 
the diffusion region (Figs. 2, 3 and 4). These regions 
do not move at a particular wave speed, but expand out 
slowly throughout the duration of the experiment. Due 
to the loss of equilibrium within the separator current 
layer, the region surrounding the current layer, which is 
non-resistive, responds in order to try and regain force 


balance by trying to rebuild the current within the sep¬ 
arator current layer. The forces present are basically 
the same as those found during the formation of the 
initial MHS equilibrium [Stevenson et a/., 2015]. 

Outside the cusp regions, the inward directed mag¬ 
netic pressure forces once more dominate over the out¬ 
ward directed plasma pressure force, generating an in¬ 
flow towards the separator current layer in this region. 
Inside the cusp regions, outward directed magnetic ten¬ 
sion forces dominate over the inward directed plasma 
pressure forces causing an outflow. These flows are 
maintained because, as soon as the current within the 
current layer strengthens to |j| = jcrit^ diffusivity dis¬ 
sipates this current and thus prevents a static equilib¬ 
rium being formed. Instead, a system, which is close to 
a steady state, is created involving slow reconnection at 
the separator current layer (phase II of the reconnection 
process). 

4.3. Infinite-Time Collapse 

Fig. 2a shows that the current enhancements on 
the separatrix surfaces are also perturbed, but there is 
no corresponding perturbation in pressure at this time 
(Fig. 2d). Furthermore, the current along the sepa¬ 
ratrix surfaces is not decreased, as one might expect if 
these perturbations were the result of reconnection, but 
instead is increased. 

Since the separatrix surfaces are outside the diffusion 
region, in an area where the magnetic plasma is non- 
resistive, the plasma continues to behave as it did in 
the relaxation experiment discussed in detail in Steven¬ 
son et al. [2015]. In particular, we recall that our ini¬ 
tial MHS equilibrium is not a true equilibrium, but is 
in force balance everywhere, except within the current 
layers on the separatrices and along the separator (the 
formation of a true equilibrium would take an infinite 
time). On the edges of the separatrix surfaces there 
are small residual forces that are very slowly increasing 
the current within the current layer and this is what 
continues to happen in our experiment here. 

4.4. Vorticity and Velocity 

In order to understand the nature of the flows cre¬ 
ated as a result of the reconnection, we consider both 
the vorticity and the velocity at three different times. 
The final row of graphs in Figs. 2, 3 and 4 show the 
vorticity in the 2 ; = 0.4 plane perpendicular to the sep¬ 
arator, through the depth and across the width of the 
current layer, respectively. Over plotted on these graphs 
are arrows indicating the direction and size of the ve¬ 
locity in these cuts. The arrows are coloured depending 
on their position: the arrows are coloured black where 
d<0 and u;<0 and the arrows are coloured green where 
d>0 and u;>0. We have coloured the arrows this way 
so that the direction of the flow from a given side of the 
diffusion region is clearer. 
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Figs. 2g-2i show a very similar pattern to the clas¬ 
sical quadrupolar vortex scenario and stagnation flow 
found in 2D X-point reconnection regimes. The main 
difference is that instead of finding zero vorticity in the 
vicinity of the separator, an antiparallel flow is found 
associated with a clockwise (blue) rotating flow pattern 
(Fig. 2i). 

If instead of considering the z = 0.4 plane, we looked 
at the vortex pattern in the z = 0.6 plane (see movie of 
Fig. 2), then we would find a similar quadrupole vortex 
(rotated slightly due to the twisted nature of the cur¬ 
rent layer and separartix surfaces about the separator), 
but in the vicinity of the separator, an antiparallel flow 
associated with an anticlockwise (red) rotation would 
be found. This agrees with the existence of a counter¬ 
rotating flow along the separator discussed in Stevenson 
and Parnell [2015]. 

Not surprisingly, therefore, looking at the vorticity in 
the cut through the depth (Fig. 3) and across the width 
(Fig. 4) of the current layer, we see that the directions 
of the flows change with position along the separator. 
The velocity arrows indicate that the dominant flows 
are directed inwards through the depth, and outwards 
across the width, for all z. However, superimposed on 
these are weak stagnation-type flow patterns. 

Along the separator, in the cut through the depth 
(Figs. 3g-3i), weak flows run towards both nulls from 
a point 0.6 times the length of the separator (as mea¬ 
sured from the lower null). The location of this stag¬ 
nation point, which corresponds to where the flows are 
purely directed inwards to the current layer through the 
depth, does not appear to move over time. Its location 
is a result of the fact that the plasma pressure on the 
separator is greatest at this point in the MHS equilib¬ 
rium and so the strongest magnetic pressure force must 
have existed at this location to counter the largest pres¬ 
sure force that would have been located there. 

The cuts across the width of the current layer 
(Figs. 4g-4i) show the opposite quadrupolar-vortex pat¬ 
tern close to the separator. This pattern shows that 
there are weak flows that run in from the nulls along 
the separator to a point 0.4 times the length of the 
separator (as measured from the lower null) at time 
t = O.OlOtf (Fig. 4g). This is the approximate centre 
of the main stagnation outflow across the width of the 
current layer. It also corresponds to the point where 
the current is largest along the separator in the MHS 
equilibrium, and thus where the outwardly directed ten¬ 
sion force must have been highest. So, as soon as the 
equilibrium is lost, a strong outflow in the z = 0.4 plane 
results. 

Unlike the stagnation flow through the depth of 
the current layer, the stagnation-point flow across the 
width changes over time, such that at the time phase I 
ends, there seem to be two stagnation points on the 
separator. During phase H a single stagnation-type 
flow has reformed. It is not clear what the multiple 


stagnation-point flows indicate, but it is important to 
remember that at the transition between the two phases 
very little reconnection, if any, is occurring since the 
slow, impulsive-bursty reconnection phase has not yet 
formed. 

5. Transport of Energy 

We have already seen that the reconnection at the 
separator current layer leads to waves being launched 
in the system, from the edges of the diffusion site, due to 
the sudden lack of force balance. These waves travel out 
and cause the magnetic field and plasma to change, set¬ 
ting up flows in the system. In this section, we analyse 
the transport of magnetic, internal and kinetic energy 
equations (as detailed in Birn et al. [2009]), integrated 
over volumes within our domain, in order to see what 
quantities these waves and flows carry with them. 

To determine the quantities in the transport equa¬ 
tions, we integrate each of them separately over sub¬ 
volumes that increase in size within our domain. The 
nine volumes have sizes —(0.15 + k/10) < x^y < 
(0.15 + k/10), for k = 0,1,2,..., 8 and the z-range is 
fixed for all volumes at —0.2 < z < 1.2, since the waves 
and flows travel horizontally and not vertically out from 
the separator. Therefore, the smallest volume encloses 
the current layer, the second volume encloses the first 
volume and so on up to the largest volume, which is 
slightly smaller than the domain size in the x and y 
dimensions. A cartoon of these volumes is shown in 
Fig. 6a, where the volumes are coloured black, pur¬ 
ple, blue, lime, green, yellow, orange and red as they 
increase in size. These volumes are shown drawn with 
the MHS equilibrium skeleton in Fig. 6b to highlight 
the size of the boxes compared to the skeleton. Hence, 
in each plot of Fig. 7, which shows the time evolution 
of the energy transport quantities, there are nine lines 
coloured to match these volumes. 

5.1. Transport of Magnetic Energy 

The transport of magnetic energy equation states 
that 

where t is time, is the square of the magnitude of the 
magnetic field {B = |B|), /tq is the magnetic permeabil¬ 
ity which is equal to one in our dimensionless units, p 
is the square of the current magnitude (j = |j|) and E 
is the electric field. Hence, the rate of change of mag¬ 
netic energy, throughout the reconnection experiment 
(Fig. 7a), is made up of the negative sum of the Ohmic 
dissipation (Fig. 7b), the Poynting flux (Fig. 7c) and 
the work done by the Lorentz force (Fig. 7d). 
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Figure 6. (a) Cartoon depicting the nine volumes over 
which the transport of energy equations are integrated. 
The volumes increase according to — (0.15 + /c/lO) < 
x,y < (0.15 + /c/lO) and —0.2 < z < 1.2 shown by the 
colours black (/c = 0), purple {k — 1), blue (k — 2)^ cyan 
(k — 3), lime (k — 4), green (k — 5), yellow (k — 6), 
orange {k = 7) and red (/c = 8), with the last box being 
just smaller than the size of the domain in x and y. (b) 
Eor context, the MHS equilibrium skeleton with boxes 
overdrawn. 

At the start of the experiment, there is an immedi¬ 
ate drop in the rate of change of magnetic energy corre¬ 
sponding to the strong Ohmic dissipation during phase 
I, the fast-reconnection phase (Fig. 7b). The integrated 
values of the Ohmic dissipation, which are of the order 
10“^, are the same regardless of the size of the volume 
(so only one line is visible) since the Ohmic dissipation 
occurs in the diffusion region, which is enclosed in all 
boxes. The Ohmic dissipation is fairly constant after 
t = 0.09t/; the slow, impulsive-bursty second phase. 

The value of the Poynting flux, integrated over all 
volumes, is initially positive indicating that the waves 
are travelling out through the boundaries of our volumes 
(Fig. 7c). The waves then cause the plasma to change, 
setting up flows in the system near to the original diffu¬ 
sion site. The flows bring Poynting flux in through the 
smaller volumes, from about t = 0.16t/ onwards. How¬ 
ever, the Poynting flux, carried out through the volumes 
by the waves, becomes relatively large the further out 
they travel (see the green to red lines in Fig. 7c). Note, 
however, that the amount of Poynting flux is roughly 25 
times smaller than the peak Ohmic dissipation in phase 
1 . 

The Lorentz force is working to try to regain force 
balance in the system from the moment the current in 
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the separator current layer begins to be dissipated (pos¬ 
itive values in Fig. 7d). This figure shows that the work 
done by the Lorentz force, which is roughly of the or¬ 
der of 3 X 10“^, is acting out through the sub-volumes 
over which we have integrated. This is directly related 
to the direction of the magnetic tension and magnetic 
pressure forces which make up the Lorentz force. The 
magnetic tension force, which acts to straighten the field 
lines, is directed outwards from the diffusion site both 
within and outwith the cusp regions which are formed 
by the separatrix surfaces of the nulls. The magnetic 
pressure force is directed in towards the diffusion region 
within the cusps and outwith the cusps. Overall, these 
forces sum such that the work done by the Lorentz force 
is acting outwards away from the diffusion region. The 
magnitude of this term is roughly 30 times smaller than 
the Ohmic dissipation term. 

5.2. Transport of Internal Energy 

Eq. 2 is the transport of internal energy equation 

^(pe) = vf - V ■ ((p + /)e)v) - (-V ■ Vp), (2) 

where p is the density and e is the internal energy per 
unit mass. 

Here, we can write pe = 3p/2, since our closure equa¬ 
tion is e = pjp{pi — 1), and 7 = 5/3 and therefore 
p + pe = bp/2. This equation states that the rate of 
change of internal energy is due to the Ohmic heating 
minus the enthalpy flux minus the work done by the 
pressure force. Fig. 7e shows the rate of change of in¬ 
ternal energy, which is of the order of 10“^, integrated 
over all nine volumes, throughout the experiment. The 
initial sharp spike in this figure is due to the Ohmic 
heating (Fig. 7b) as was seen in the energetics discussed 
in Stevenson and Parnell [2015]. 

After this spike, the rate of change of internal en¬ 
ergy decreases and becomes negative. A small travelling 
wave is seen moving out through all the sub-volumes 
which comes from the enthalpy flux term (Fig. 7f). 
This term is of the order of 5 x 10“^ and is half the size 
of the peak Ohmic dissipation during phase 1. 

The final term which contributes to Eq. 2 is the work 
done by the pressure force (negative values in Eig. 7d), 
however, the magnitude of this term is about 30 times 
smaller than that of the Ohmic heating and about 15 
times smaller than the enthalpy flux and so its contri¬ 
bution is small here. The work done by the pressure 
force is directed in through the sub-volumes over which 
we integrate and is acting, like the Lorentz force, to 
try to regain force balance in the system as soon as the 
reconnection begins. 
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Figure 7. Quantities, plotted against time, of (a) the rate of change of the magnetic energy, (b) the 
Ohmic dissipation, (c) the Poynting flux, (d) the work done by the Lorentz, viscous and pressure forces, 
(e) the rate of change of internal energy and (f) the enthalpy flux. The line colour represents the different 
volumes over which these quantities have been integrated (c.f. Eig 6). The black dashed vertical line 
highlights where phase I ends and phase II begins and the black dashed horizontal line indicates zero. 


5.3. Transport of Kinetic Energy 

The transport of kinetic energy equation states that 
the rate of change of kinetic energy is equal to the work 
done by the Lorentz force plus the work done by the 
pressure force plus the work done by the viscous force 
minus the bulk kinetic energy flux 

= v-(jxB) + (-v-Vp) + v-F^-V - 

(3) 

where = z^(V^v + |V(V • v)) is the viscous force. 

The rate of change of kinetic energy is very small 
5 X 10“^) throughout the reconnection experiment. 
This is because the work done by the Lorentz and pres¬ 
sure forces (Figs. 7d) are about equal in size, but are 
of opposite sign (they are acting to regain force balance 
in the system after the dissipation of the current layer). 
Also, contributions from the work done by the viscous 
force (values close to zero in Fig. 7d) and the bulk ki¬ 
netic energy flux are very small (^ 5 x 10“^) since the 
velocities in the system have small magnitudes. 

Overall, we have found that there are five main terms 
which play a significant role in the transport of energy 
in our experiment. Ohmic heating plays the most sig¬ 


nificant role in our experiment, especially during phase 
I, converting magnetic energy into internal energy. This 
energy is then carried away from the diffusion region by 
the enthalpy flux, which is half the size of the Ohmic 
heating term and the Poynting flux, which is roughly 25 
times smaller than the peak Ohmic heating. The final 
two important terms are the work done by the Lorentz 
and pressure forces, which are similar in magnitude, but 
act in opposite directions and are both about 30 times 
smaller than the Ohmic heating term. 

6. Conclusions 

In this paper, we have studied the properties of the 
waves and flows created due to spontaneous reconnec¬ 
tion at a 3D separator current layer. We start with a 
system that is in MHS equilibrium everywhere save for 
very small forces at the current enhancements about the 
separator and separatrix surfaces. An anomalous diffu- 
sivity is applied such that reconnection only occurs at 
the separator current layer that twists about the sepa¬ 
rator. 

The onset of the reconnection produces waves that 
propagate out from the edge of the diffusion site at 
the separator current layer. These waves only have 
small amplitudes, due to the relatively small recon- 
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nection event that initiates them, and they travel at 
the fast-magnetoacoustic speed (which in our high- 
beta experiment is approximately equal to the slow- 
magnetoacoustic speed). 

They carry the dissipated current away from the dif¬ 
fusion region and disperse as they travel. The nature of 
the waves has the same pattern in all planes perpendic¬ 
ular to the separator, which is basically the same as that 
found due to waves launched following reconnection at 
a 2D X-type null. 

(i) Planar-like waves that are twisted about the 
separator are launched from either side of the diffu¬ 
sion region and travel away from the separator current 
layer carrying current and causing a deficit in pressure. 
Equivalent waves are also launched inwards through the 
depth of the current layer at the same time. These 
waves end up running closely behind the outwardly- 
launched waves. 

(ii) Point-like waves that again are twisted about the 
separator are launched outwards from the narrow edges 
of the diffusion region. In any given plane perpendicu¬ 
lar to the separator, these spread in a circular pattern 
carrying current away from the separator current layer 
and causing an enhancement in pressure. As above, 
point-like waves also travel inwards across the width of 
the separator current layer, which is roughly 20 times 
the size of the depth, thus these waves lag behind the 
outward waves. 

These waves communicate the sudden loss of force 
balance within the separator current layer and, hence, 
in their wake magnetic and plasma forces are set up 
with the aim of restoring the equilibrium. As already 
explained by Stevenson et al [2015], an equilibrium in 
such a system with a separator involves a current en¬ 
hancement about the separator and, thus, a velocity 
flow pattern is created, which brings in more flux from 
outwith the cusp regions to enhance the current at the 
separator. As soon as the current in the layer reaches 
the level of jcrit the anomalous diffusivity dissipates it. 
This leads to a slow, impulsive-bursty second phase of 
reconnection driven by this stagnation-like flow. Al¬ 
though the reconnection during this phase is very slow, 
the Ohmic dissipation associated with it is still larger 
than the viscous heating associated with the damping 
of the magnetoacoustic waves and flows. 

The amplitude of the waves that result from the re¬ 
connection are small and of the order 10“^. They are, 
however, much bigger (100 times) than those found 
in the 2D high-beta experiment of Fuentes-Fernandez 
et al [2012a] which we believe is due to the third di¬ 
mension permitting a larger current layer to be formed. 
To form an even larger current layer in our high-beta 
scenario we could have start, prior to forming the MHS 
equilibrium, with an initial magnetic field that has ei¬ 
ther greater initial current or different magnetic field 
parameters (see Stevenson, J. E. H. [2015] for full de¬ 
tails). However, in the resulting MHS equilibrium the 


current layer was not resolved and we were concerned 
that numerical diffusion had occurred, in some cases, 
before the MHS equilibrium could be formed. Also, the 
resulting waves were only marginally greater in ampli¬ 
tude. Lowering the value of jcrit ^ the level above which 
the diffusivity is non-zero, would define a larger cur¬ 
rent layer, but this has the side effect of permitting re¬ 
connection within the current enhancements along the 
separatrix surfaces leading to more (complicated) wave 
pulses. 

An analysis of the energy transport in the model 
shows that the Ohmic dissipation is about twice that of 
the enthalpy flux carried by the magnetoacoustic waves 
away from the reconnection site. In turn, the enthalpy 
flux is more than ten times the work done by either the 
Lorentz or the pressure forces (which are basically the 
same size, but cancel each other out since they work in 
opposite directions). The Poynting flux is also about 
ten times smaller than the enthalpy flux. The dom¬ 
inance of the enthalpy flux over the Poynting flux is 
not surprising since our experiments are high beta (as 
shown in Birn et al [2009]). 

In order to compare our dimensionless results to those 
of a specific space plasma scenario, the speed of the 
waves generated by the reconnection can be scaled by 
the factor B ^/In situ measurements by Double 
Star have detected number densities of = 10 x 10^ 
m“^ and magnetic field strengths of Bn = 30 nT at 
the magnetopause [Trenehi et al, 2008]. Using a mean 
particle mass (this value has been calculated using mag¬ 
netosheath abundances [Gloeekler and Hamilton, 1987] 
of the major magnetospheric ions [Reme et al, 2001]) 
of m = 1.07 X rup, where rup is the mass of a proton, 
this number density corresponds to a mass density of 

= 1.8 X 10“^^ kgm“^. In our experiment, the maxi¬ 
mum dimensionless Alfven speed in the outflow regions 
is va = 1-6, thus we find that the maximum Alfven 
speed in our domain is va = 320 kms“^. This is of 
the same order as the hybrid Alfven speed, vah = 380 
kms“^, found between the magnetosheath and the mag¬ 
netosphere in Komar et al. [2013]. Note that, in our 
model, the density in the outflow region decreases away 
from the separator and the magnetic field strength in¬ 
creases, so if our domain was larger we would be able 
to measure greater values of the Alfven speed than we 
have detailed here. 

All the experiments discussed here can only be run 
until the travelling waves near the boundaries of the 
box. If it was possible to run these experiments for 
longer then the viscous heating may increase sufficiently 
to become comparable with the observed Ohmic heating 
in this second phase. Furthermore, a low-beta system 
may also permit greater viscous heating. We plan, in a 
follow up paper, to investigate if this is the case. 
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